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Abstract 

We study the overlap and the fixed point Dirac operators for mas- 
sive fermions in the two-flavor lattice Schwinger model. The masses of 
the triplet (pion) and singlet (eta) bound states are determined down 
to small fermion masses and the mass dependence is compared with 
various continuum model approximations. Near the chiral limit, at 
very small fermion masses the fixed point operator has stability prob- 
lems, which in this study are dominated by finite size effects, however. 
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1 Motivation and Introduction 



Lattice Dirac operators obeying the Ginsparg- Wilson condition £Q in its sim- 
plest form 

75^ + ^75 = aDj 5 D, (1) 

have the advantage to violate chiral symmetry only locally, with quasi au- 
tomatic 0(a) improvement. Their eigenvalues lie on a circle in the complex 
plane; zero eigenvalues correspond to chiral eigenstates and indicate topolog- 
ical (instanton) modes in the gauge field. We known of at least two explicit 
incarnations of such actions, the so-called overlap action 2j and the perfect 
action Both are technically difficult and expensive to determine and to 
incorporate in a lattice simulation with dynamical fermions. 

The Schwinger model (2D fermion-gauge theory with U(l) gauge group) 
generalized to e.g. 2 flavors provides an attractive model to study such ac- 
tions. It has similar features like QCD with 2 quark flavors - although due 
to a quite different underlying dynamics. It has only bosonic bound states, 
including a massless triplet, and a non- vanishing fermion condensate. The 
overlap action may be constructed from an arbitrary action, like the Wilson 
action. It involves a computation of the sign function for hermitian matrices, 
which has to be repeated for each gauge field and is computationally chal- 
lenging. A fixed point action (the perfect action determined at the fixed point 
of the continuum limit) on the other hand has a large number of couplings 
and is difficult to construct. For the Schwinger model such a fixed point ac- 
tion has been determined in excellent approximation a few years ago jl], at 
at time, when the spectacular spectral properties of GW-fermion operators 
were not yet widely discussed. Both, the overlap and the fixed point action 
then have been studied in the Schwinger model for various aspects, mainly 
related to the Dirac operator spectrum [H] and to the topological content of 
the gauge fields. 

Introducing mass in the Schwinger model is a subtle problem [U]-|15j 
and the continuum approximations involve plausible but non-stringent ad- 
hoc assumptions Eij. Of particular interest is the analog of the PCAC- 
relation, i.e. the dependence of the triplet bound state mass on the fermion 
mass parameter. Recent studies for staggered action [Tfi] . Wilson action 
[TTj and overlap action JH] have confirmed the exponent 2/3. Here we try 
to obtain results closer to the chiral limit for two GW-type actions. The 
motivation is on one hand further confirmation of the functional dependence 
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and a clearer determination of the correct coefficient. 

On the other hand we want to compare the relative merits of these chirally 
improved actions closer to the chiral limit, which is hardly accessible with e.g. 
Wilson-type actions. In fact, actions where (near the origin) the eigenvalue 
spectra scatter have problems at small fermion masses in general. Due to 
this fluctuation of the small real eigenvalues many configurations will mimic 
zero modes and necessitate special techniques to deal with them. The hope 
is to improve this situation with overlap or perfect actions. We discuss our 
conclusion based on this model study at the end. 

2 The massive 2-flavor Schwinger model 

We study the euclidean 2d gauge theory with gauge group U(l) and two 
flavors of fermions with degenerate mass. For vanishing fermion masses, the 
classical theory has a symmetry group U(2)lxU (2)r that is broken down by 
the anomaly to SU(2)l x SU {2)r x U (l)y. For non-vanishing fermion masses 
the chiral symmetry is broken at the classical level and a U(2)y symmetry 
remains. 

Let us briefly discuss some known analytical results for the continuum 
model. The generalized Sine-Gordon model bosonizes the 2-flavor Schwinger 
mo del [I]. Its Lagrangian reads 

- (V 0) ) 2 + - (V 3) ) 2 + - — (v? (0) ) 2 -— cos (V2^ (0) ) cos (V2^ (3) ) • 

2 '2 2 7T 7T 

(2) 

The constant c is not determined by the bosonization but is related (see 
e.g. [T2*l IT4*| ) to the masses /i^ and /i^ 3 -* used for normal ordering of the 
fields <^ (0) and y? (3) by c = ( / u (0) / u (3) ) 1/2 e 7 /2, where 7 (= 0.577216 . . .) denotes 
Euler's constant. For the case of two flavors two of the currents - which we 
measure - (compare (|14jl ) are bosonized with the following prescription 

r a <g> 7m V0*0 = -^e^d v v (a) (x) , a = 0,3, (3) 

where the generators r a of rotations in flavor space are given by the Pauli 
matrices r a = a a , a = 1,2,3 and r° = 1. For the other two members 
j^ij^f 1 of the iso-triplet no explicit abelian bosonization is known; however 
due to invariance under flavor rotations, their masses are the same as for the 
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triplet current (a = 3, which we called pions due to the analog to QCD). For 
vanishing quark mass m in (J2J) the two flavor Schwinger model is bosonized 
by two free fields. One of them, the pion field ^ 3 \ is massless (M n = 0) 



and the eta field (p^ obtains the Schwinger mass M v = gy2/-K (due to the 
anomaly) . 

For non-vanishing quark mass, also the bosonized model (J2J) can no longer 
be solved in closed form. A semi-classical analysis (see e.g. [Hj) should pro- 
vide a good approximation when all involved masses are large. The squares 
of the masses are given by the second derivatives of the interaction part 
^ (0) ,<^ (3) ] of © at the minimum. After normal ordering the fields <p^ a > 
with respect to their own masses (setting //°) = M v and ^ = M n ), the 
semi-classical analysis for the iso-triplet and iso-singlet gives 



Mr, 



\ 




(5) 



The first relation plays a role similar to the PCAC relation in QCD (where 
the pion mass squared is proportional to the quark mass). 

In an attempt to go beyond semi-classical analysis one can approximate 
the generalized Sine-Gordon model Q by a solvable model. In the limit of 
large coupling [7] g and small mass m the iso-singlet field becomes static 
and the model is reduced to a standard Sine-Gordon model for the triplet 
field (p^ 3 \ with the Lagrangian 

i(V 3) ) 2 - 2Ccos(v / 2^ {3) ) . (6) 

This reduced bosonized theory has been studied HHj using the WKB ap- 
proximation [2D] and recently again by Smilga ^21- His analysis is now based 
on the newly derived analytic solution [21] of the standard Sine-Gordon model 
and he finds 

(7) 

The numerical coefficient of the (m/g) 2 ^ 3 term is 2.008, somewhat smaller 
than the value 2.1633 in (jlj). The truncation of the original model (j2J to 
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the standard form (jHJ) does not allow for a result for the singlet mass M v . 
For this state only the semi-classical formula (JSJ) is available. It his however 
interesting to use the exact result (J7|) of the truncated model as an input in 
(jSj), and below we compare also this formula to the numerical data. 

3 Lattice Dirac Operators 

The Schwinger model is a super-renormalizable theory, i.e. the bare coupling 
does not get renormalized and equals the physical gauge coupling g. On 
the lattice for the standard Wilson plaquette gauge action we use the usual 
coupling f3 in front of the gauge field action. In the naive continuum limit 
this is related to the continuum coupling via (3 = 1/g 2 (we set the lattice 
spacing a = 1). We define the coupling g which we use for comparison with 
the analytical results in the continuum (see above) through this relation. 

We first generate the massless fermion actions and introduce the fermion 
mass parameter subsequently. 

Overlap Dirac operator 

For the overlap action [22] one may start with the Wilson Dirac operator 

Avi (x, y) = {m w +2) 5 xy -- ^ [(1 + 07J U xy 5 x , y -^ + (1 - c^) U\ x 6 XfV +p , 

(8) 

(with the Pauli matrices a^) at some value of m w G (—1,0) and then con- 
struct 

D 0v = - [l + 75 e(7 5J D Wi ) ]. (9) 

Some words about the choice of m w : according to (221 ^ ^ s arbitrary, in the 
sense that any (strictly negative) value in the interval (—1,0) reproduces 
the correct continuum theory, but it may be optimized with regard to its 
scale dependence by looking for example at the behavior of the (projected) 
spectrum. In many practical determination of the sign function configura- 
tions with small eigenvalues of 75 D Wi are cumbersome; this situation may 
be improved by a suitable choice of m w . Comparing expectation values of 
operators like (ipi/j) for different m w one has to take care of the proper nor- 
malization 23J. Comparing with free lattice fermions one finds a (trivial) 
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factor of \J1m w renormalizing the field operator, i.e. (^ip) = Zl {ipip) ov 
with Zj/j = \2m w \ in our convention. We choose m w — — 1. 

The operative definition of the generalized sign function e( 7 5 Dyj\ ) enter- 
ing the above equation is through its eigenvalues, 

e ( 75 D m ) = U Sign(A) C/t with 7 5 D Wi = U Arf . (10) 

Here Sign(A) denotes the diagonal matrix containing the signs of the eigenval- 
ues obtained through the unitary transformation U of the hermitian matrix 
75 Dwi • There are various ways to numerically find D<j v without passing 
through the diagonalization problem (which is prohibitively expensive for 
D — 4) [2H 123 I2H] (for D = 2 see also [21] and recently, in the Schwinger 
model [23 HE])- I n our simple context computer time is no real obstacle and 
therefore we use the direct definition ([10)1. explicitly performing the diagonal- 
ization. For comparison we also applied an iterative technique implementing 
Newton's method [2UI2ZI- For many gauge configurations the resulting oper- 
ators agreed to the requested accuracy (8 digits). For larger (3 > 4 however, 
the number of configurations with convergence problems for the iterative 
method due to small eigenvalues increased. The results presented here are 
all based on the exact diagonalization. 

Fixed point Dirac operator 

In Ref.jlj the fixed point Dirac operator was parameterized as 

^(^)4eem/)^^/)> withy +5/. (ii) 

Here / denotes a closed loop through x or a path from the lattice site x 
to y = x + Sf (distance vector 5f) and U(x, f) is the parallel transporter 
along this path. The <Tj-matrices denote the Pauli matrices for % = 1, 2, 3 and 
the unit matrix for i = 0. Note the factor 1/2 in ()11|) in order to obtain the 
same normalization as the overlap Dirac operator. The action obeys the usual 
symmetries as discussed in j3]; altogether it has 429 terms per site. The action 
was originally determined at large (3 for gauge fields distributed according to 
the non-compact formulation with the Gaussian measure. Excellent scaling 
properties, rotational invariance and continuum-like dispersion relations were 
observed at various smaller values of the gauge coupling (3. 
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In Ref . |Hj this action was studied both, for compact and the original non- 
compact gauge field distributions. In the compact case the action is not 
expected to exactly reproduce the fixed point of the corresponding BST, 
but nevertheless it is still a solution of the GWC violations are instead 
introduced by the parameterization procedure, which cuts off the less local 
couplings. It was demonstrated, that the eigenvalue spectrum is close to 
circular shape, improving in this respect towards larger (5 (see [2E] for a more 
detailed discussion). 

Here we study the action only for the compact gauge field distributions 
in order to allow for a direct comparison with the overlap Dirac operator. 

The massive case 

The massive overlap Dirac operator may be related to the massless one by 



where the mass parameter — 1 < /j, < 1 is related to the fermion mass by 



with Z^Z m = 1, where Z m and Z$ are the mass and the wavefunction renor- 
malization constants, respectively. Since with our choice of parameterization 
the massless fixed point Dirac operator (fTTj) has the same overall renormal- 
ization as the massless overlap Dirac operator we also use ffT^j) to generate 
the massive fixed point Dirac operator. 

4 Simulation Details 

Uncorrelated gauge configurations for lattice size 16 2 and 24 2 have been gen- 
erated in the quenched setup. However, we are including the fermionic de- 
terminant in the observables: all the results presented here are obtained with 
the correct determinant (squared, for two flavors) weight. From earlier expe- 
rience [HE] we learned that this is justifiable for this model and the presented 
statistics. We perform our investigation on two sets of 5000-10000 configu- 
rations at /3 — 4 and 6. These are the same gauge configurations as used in 
Ref . [T7j for the determination of bound state masses for the Wilson Dirac 
operator and thus we may compare the relative efficiency in the results. 
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The configurations have been well separated by three times the auto- 
correlation length for the so-called "geometric definition" of the topological 
charge. For each configuration we construct D Qv and D Fp as discussed. For 
each lattice Dirac matrix we then compute the inverse (the quark propagator) 
and the determinant. 

The masses for the pions and the eta-particle were determined from the 
decay of two point functions of the following vector currents 

J / = ^ x 1 ® 7/1 V'* for the//. (14) 

We emphasize that 2-point functions of scalars or pseudo-scalars are not 
particularly well suited for the determination of the meson masses. These 
operators are bosonized (compare the discussion above) by cosines and sines 
of the fundamental fields ^] and their 2-point functions strongly mix 
contributions from both the triplet and singlet states. For the two point 
functions of the momentum-zero projected vector currents one has to take 
care about the proper choice of \i. As could be verified by gauge field inte- 
gration there are no contributions to the propagation of J% 2 (p = 0) in the 
direction \i = 2. 

5 Results 

Let us now discuss our numerical results for the mass spectrum. We present 
the data for the overlap Dirac operator for j3 = 4 and 6 and L = 16 and 24 
in Fig. ^ The results for the fixed point operator at /3 = 6 and L = 16, 24 
are shown in Fig. |2j In each plot we combine the values of the singlet mass 
and the triplet mass and compare the numerical results to Smilga's formula 
(JZj) (full lines) and the semi-classical results (jlj) and (jHJ (dotted lines). 

We observe that the data are well described by the semi-classical result 
(@J) (dotted line). In particular at quark masses above (m^) 2 / 3 > 0.3 the 
semi-classical approximation is fulfilled best. In this region our results for 
the triplet mass coincide with the results of [IB]. For small quark masses 
we observe a strong and seemingly systematic deviation especially for small 
lattice size. We attribute this to finite size effects. In particular, the bo- 
son propagator does not really become asymptotic at such large correlation 
lengths. 
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Figure 1: Results for overlap Dirac operator {3 — 4, 6, L — 16, 24. Symbols: 
Monte Carlo data. Full line: Smilga's formula ((Zj). Dashed line: the semi- 
classical approximation (jlj) and (0). 
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Figure 2: Results for the fixed point Dirac operator at j3 = 6, L = 16, 24. 
Symbols: Monte Carlo data. Full line: Smilga's formula (jJJ). Dashed line: 
the semi-classical approximation (J3J) and (jSJ. 



For the fixed point operator we only get reasonably results for small m 
at large gauge couplings j3 > 4, when the spectrum is close to circular. For 
smaller /3 the scattering of the small eigenvalues is too large and of the order 
of the bare fermion mass parameter, causing disturbing fluctuations. The 
dispersion a of the eigenvalues around the circle (which decreases with larger 
(3[29 J leads to accidental poles in the propagators even at non- vanishing 
quark mass. 

For the eta mass M n we find that the data approach the correct m = 
value as the quark mass vanishes. The semi-classical formula (jSJ provides a 
reasonable description of the data. For small m the quantum fluctuations 
become larger and thus the numerical data deviate from the semi-classical 
curve. The cases with comparatively large error bars always indicate lack of 
sufficient statistics for the propagator. The behavior at small m can be fitted 
by a power law as has been done in Ref . [TK] for the results from the model 
with staggered fermions. 

Comparing the results with those obtained in Ref.[T7] for the Wilson 
Dirac operator we find generally smoother behavior, in particular at smaller 
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fermion mass parameters. This was to be expected due to the improved 
chirality properties. On the other hand, for a full exploitation of this feature 
larger lattices with better finite size control are needed. 

6 Conclusion 

We find that our results for the singlet mass and the triplet mass are well 
described by the semi-classical formulas (Jlj). (j3j). For the triplet mass M n 
we observe a clear deviation from the semi-classical curve at small quark 
masses which we attribute to finite volume effects. For the 77-mass the data 
are in good agreement with the semi-classical approximantion with quantum 
fluctuations becoming less important as the quarks are made heavier. 

As for the studied model, although the results for 77 and 7r are in good 
agreement with the the analytic expectations the values of the pion mass 
below < 0.3 are clearly plagued by the limited lattice size. However, 

our emphasis was on the study of the effects due to the chiral properties of 
the Dirac operators. For better results for the 7r-mass at very small masses 
much larger lattice and higher statistics appear to be necessary to allow a 
control of the finite size effects. 

For small fermion masses it seems to be important, that the eigenvalues 
of the Dirac operator do not scatter too much away from the optimal unit 
circle. With this regard the overlap action does the better job. Even the 
qualitatively quite good fixed point action exhibits a weakness for small (3 due 
to this effect. Spurious singularities or almost-singularities in the propagator 
necessitate very large samples in order to obtain reliable mean expectation 
values and thus propagators. The width of the distribution of eigenvalues 
around zero therefore essentially limits the practically applicable values of 
the bare fermion mass parameter. 

Acknowledgment: We want to thank C. Gattringer and I. Hip for many 
useful discussions. 
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